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Anderson Orthogonality (AO) refers to the fact that the ground states of two Fermi seas that 
experience different local scattering potentials, say Id) and |Gf), become orthogonal in the thermo- 
dynamic limit of large particle number N, in that |(Gi|Gp)| ~ ]\[~2'^ao for N — >■ oo. We show that 
the numerical renormalization group offers a simple and precise way to calculate the exponent Aaq: 
the overlap, calculated as function of Wilson chain length k, decays exponentially, ~ e~*", and 
Aaq can be extracted directly from the exponent a. The results for Aaq so obtained are consistent 
(with relative errors typically smaller than 1%) with two other related quantities that compare how 
ground state properties change upon switching from |Gi) to |Gp); the difference in scattering phase 
shifts at the Fermi energy, and the displaced charge flowing in from infinity. We illustrate this for 
several nontrivial interacting models, including systems that exhibit population switching. 

PACS numbers: 



I. INTRODUCTION 

In 1967, Anderson considered the response of a Fermi 
sea to a change in local scattering potential and made the 
following observation:^ the ground states |Gi) and |Gf) of 
the Hamiltonians Hi and Hp describing the system before 
and after the change, respectively, become orthogonal in 
the thermodynamic limit, decaying with total particle 
number N as 

|(Gi|GF)|^iV-5Aio , (1) 

because the single-particle states comprising the two 
Fermi seas are characterized by different phase shifts. 

Whenever the Anderson orthogonality (AO) exponent 
Aao is finite, the overlap of the two ground state wave 
functions goes to zero as the system size becomes macro- 
scopic. As a consequence, matrix elements of the form 
KGiIOIGf)!, where O is a local operator acting at the 
site of the localized potential, necessarily also vanish in 
the thermodynamic limit. This fact has far-reaching con- 
sequences, underlying several fundamental phenomena 
in condensed matter physics involving quantum impu- 
rity models, i.e. models describing a Fermi sea coupled 
to localized quantum degrees of freedom. Examples are 
the Mahan exciton (ME) and the Fermi-edge singularity 
(FES)'^"'' in absorption spectra, and the Kondo effect*^' 
arising in magnetic alloys ' or in transport through quan- 
tum dots.^ For all of these, the low-temperature dynam- 
ics is governed by the response of the Fermi sea to a 
sudden switch of a local scattering potential. More re- 
cently, there has also been growing interest in inducing 
such a sudden switch, or quantum quench, by optical 
excitations of a quantum dot tunnel-coupled to a Fermi 
sea, in which case the post-quench dynamics leaves fin- 
gerprints, characteristic of AO, in the optical absorption 
or emission line shape. '^^^^ 

The intrinsic connection of local quantum quenches to 
the scaling of the Anderson orthogonality with system 
size can be intuitively understood as follows. Consider 



an instantaneous event at the location of the impurity at 
time t = in a system initially in equilibrium. This local 
perturbation will spread out spatially, such that for t > 0, 
the initial wave function is affected only within a radius 
L V ft of the impurity, with Vf the Fermi velocity. The 
AO finite-size scaling in Eq. (1) therefore directly resem- 
bles the actual experimental situation, and in particular 
allows the exponent Aao to be directly related to the 
exponents seen in experimental observables at long time 
scales, or at the threshold frequency in Fourier space. 

A powerful numerical tool for studying quantum im- 
purity models is the numerical renormalization group 
(NRG),^'''^'* which allows numerous static and dynami- 
cal quantities to be calculated explicitly, also in the ther- 
modynamic limit of infinite bath size. The purpose of 
the present paper is to point out that NRG also offers a 
completely straightforward way to calculate the overlap 
|(Gi|Gf)| and hence to extract Aao- The advantage of 
using NRG for this purpose is that NRG is able to deal 
with quantum impurity models that in general also in- 
volve local interactions, which are usually not tractable 
analytically. Although Anderson himself did not include 
local interactions in his considerations,^ his prediction (1) 
still applies, provided the ground states |Gi.f) describe 
Fermi liquids. This is the case for most impurity models 
(but not all; the two-channel Kondo model is a notable 
exception). Another useful feature of NRG is that it al- 
lows consistency checks on its results for overlap decays, 
since Aao is known to be related to a change of scat- 
tering phase shifts at the Fermi surface. These phase 
shifts can be calculated independently, either from NRG 
energy fiow diagrams, or via Friedel's sum rule from the 
displaced charge, as will be elaborated below. 

A further concrete motivation for the present study is 
to develop a convenient tool for calculating AO exponents 
for quantum dot models that display the phenomenon of 
population switching. "'^"■^■^ In such models, a quantum 
dot tunnel-coupled to leads contains levels of different 
widths, and is capacitively coupled to a gate voltage that 
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shifts the levels energy relative to the Fermi level of the 
leads. Under suitable conditions, an (adiabatic) sweep of 
the gate voltage induces an inversion in the population 
of these levels (a so called population switch), implying 
a change in the local potential seen by the Fermi seas 
in the leads. In this paper, we verify that the method 
of extracting Aao from (GilGp) works reliably also for 
such models. In a separate publication^^ we will use this 
method to analyze whether AO can lead to a quantum 
phase transition in such models, as suggested in Ref. 23. 

The remainder of this paper is structured as follows: In 
Sec. II we define the AO exponent Aao in general terms, 
and explain in Sec. Ill how NRG can be used to calculate 
it. Sec. IV presents numerical results for several interact- 
ing quantum dot models of increasing complexity: first 
the spinless interacting resonant level model (IRLM), 
then the single- impurity Anderson model (SIAM), fol- 
lowed by two models exhibiting population switching, 
one for spinless, the other for spinful electrons. In all 
cases, our results for Aao satisfy all consistency checks 
to within less than 1%. 



II. DEFINITION OF ANDERSON 
ORTHOGONALITY 

A. AO for a single channel 

To set the stage, let us review AO in the context of 
a free Fermi sea involving a single species or channel of 
noninteracting electrons experiencing two different local 
scattering potentials. The initial and final systems are 
described in full by the Hamiltonians Hi and Hp, respec- 
tively. Let c^xlO) be the single-particle eigenstates of 
Hx characterized by the scattering phase shifts Sx{£), 
where X G {I, F} and are fermion creation opera- 
tors, and let e-^ be the same Fermi energy for both Fermi 
seas \Gx)- Anderson showed that in the thermodynamic 
limit of large particle number, N — )■ oo, the overlap 

(Gi|Gp) = (0| n n <fIO) (2) 

e<ef e<ef 

decays as in Eq. (1),^'"' where Aao is equal to the differ- 
ence in single-particle phase shifts at the Fermi level, 

AAo = Aph = [(5F(eO-'5i(£0]A- (3) 

The relative sign between Aao and Aph (+, not — ) 
does not affect the orthogonality exponent A^q , but fol- 
lows standard convention [Ref. 24, Eq. (7), or Ref. 25, 
Eq. (21)]. 

In this paper we will compare three independent ways 
of calculating Aao- (i) The first approach calculates the 
overlap |(Gi|Gf)| of Eq. (1) explicitly as a function of 
(effective) system size. The main novelty of the present 
paper is to point out that this can easily be done in the 
framework of NRG, as will be explained in detail in sec- 
tion III. 



(ii) The second approach is to directly calculate Aph 
via Eq. (3), since the extraction of phase shifts Sx{£-^) 
from NRG finite-size spectra is well-known: ^'^ provided 
that Hx describes a Fermi liquid, the (suitably normal- 
ized) fixed point spectrum of NRG can be reconstructed 
in terms of equidistant free-particle levels shifted by an 
amount determined by Sx{s^). The many-body exci- 
tation energy of an additional particle, a hole and a 
particle-hole pair thus allow the phase shift 6xi£-^) to 
be determined unambiguously. 

(iii) The third approach exploits Friedel's sum rule,^^ 
which relates the difference in phase shifts to the so called 
displaced charge Ach via Ach = Aph- Here the displaced 
charge Ach is defined as the charge in units of e {i.e. the 
number of electrons) flowing inward from infinity into a 
region of large but finite volume, say Marge, surrounding 
the scattering location, upon switching from Hi to Hp: 

Ach = (GF|ntot|GF)-(Gi|ntot|Gi) 

= A,ca + Adot • (4) 

Here ntot = 'T^sca + ^^-dot , where ngoa is the total number of 
Fermi see electrons within t^argc, whereas Udot is the local 
charge of the scattering site, henceforth called "dot" . 
To summarize, we have the equalities 

Aio = All, = Ac\ , (5) 

where all three quantities can be calculated indepen- 
dently and straightforwardly within the NRG. Thus 
Eq. (5) constitutes a strong consistency check. We will 
demonstrate below that NRG results satisfy this check 
with good accuracy (deviations are typically below 1%). 

B. AO for multiple channels 

We will also consider models involving several indepen- 
dent and conserved channels {e.g. spin in spin-conserving 
models). In the absence of interactions, the overall 
ground state wave function is the product of those of 
the individual channels. With respect to AO, this triv- 
ially implies that each channel adds independently to the 
AO exponent in Eq. (1), 

Aio Aio,,, (6) 

where fj, = 1, . . . ,Nc labels the A'c different channels. 
We will demonstrate below that the additive character in 
Eq. (6) generalizes to systems with local interactions, pro- 
vided that the particle number in each channel remains 
conserved. This is remarkable, since interactions may 
cause the ground state wave function to involve entan- 
glement between local and Fermi sea degrees of freedom 
from different channels. However, our results imply that 
the asymptotic tails of the ground state wave function far 
from the dot still factorize into a product of factors from 
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individual channels. In particular, we will calculate the 
displaced charge for each individual channel, cf. Eq. (4), 

Ach,Ai = (Gf| n-tot,M I^'f) - (Gil ntot,p |Gi) 

= Asca,^ + Adot,^ , (7) 

where ntot,ii = j^sca,/j + j^-dot./^- Assuming no interactions 
in the respective Fermi seas, it follows from Friedel's sum 
rule that A\q ^ = A^j^ ^, and therefore 

AT, 

Aio = E^ch.,-Ac\, (8) 

where A^j^ is the total sum of the squares of the displaced 
charges of the separate channels. Equation (8) holds with 
great numerical accuracy, too, as will be shown below. 

III. TREATING ANDERSON 
ORTHOGONALITY USING NRG 

A. General impurity models 

The problem of a noninteracting Fermi sea in the pres- 
ence of a local scatterer belongs to the general class of 
quantum impurity models treatable by Wilson's NRG.^'* 
Our proposed approach for calculating Aao applies to 
any impurity model treatable by NRG. To be specific, 
however, we will focus here on generalized Anderson im- 
purity type models. They describe Nc different (and 
conserved) species or channels of fermions that hybridize 
with local degrees of freedom at the dot, while all inter- 
action terms are local. 

We take both the initial and final {X £ {I, F}) Hamil- 
tonians to have the generic form Hx = H]j + x +^int- 
The first term, 

^b = ^^ect^c,^, (9) 

/A— 1 e 

describes a noninteracting Fermi sea involving chan- 
nels. (A^c includes the spin index, if present.) For sim- 
plicity, we assume a constant density of states p^ie) = 
po^f^9{D — \e\) for each channel with half-bandwidth D. 
Moreover, when representing numerical results, energies 
will be measured in units of half-bandwidth, hence D := 
1. The Fermi sea is assumed to couple to the dot only 
via the local operators /q^^ = ^^7= and /g^, that, 

respectively, annihilate or create a Fermi sea electron of 
channel fi at the position of the dot, r = 0, with a proper 
normalization constant A'b to ensure [/o/i,/o^/] = ^tifj.'- 

The second term, H^ x, contains the non-interacting 
local part of the Hamiltonian, including the dot-lead hy- 
bridization, 

Ha.x - E ^d^.xfid^, + E V [^VoM + h-^] (-10) 

/i— 1 ju— 1 



Here ed^,x is the energy of dot level (jl in the initial or final 
configuration, and n^^ = d^^dfj, is its electron number. 
Fp = TTPfj.V^ is the effective width of level p induced by 
its hybridization with channel p of the Fermi sea, with 
the /i-conserving matrix element connecting the d-level 
with the bath states c^^, taken independent of energy, 
for simplicity. 

Finally, the interacting third term is given in the case 
of the single-impurity Anderson model (SIAM) by the 
uniform Coulomb interaction U at the impurity, 

i?™ = i[/nd(nd-l), (11) 

with rid = Xl^i^d/i, while in case of the interacting res- 
onant level model (IRLM), the interacting part is given 

by 

A^^L^ = U'n^no , (12) 

with no = E^/o.m/o.m = E^"o,ai- In particular, most 
of our results are for the one- or two-lead versions of the 
SIAM for spinful or spinless electrons, 

AM ^ ^ ^ ^ ^SIAM (13) 

We consider either a single dot-level coupled to a single 
lead (spinfuU, = 2 : p G {f, i}), or a dot with two 
levels coupled separately to two leads (spinless, Nc — 2 : 
p e {1,2}; spinfuU, A^c = 4 : p € {If, U, 2t, 2^}). A 
splitting of the energies Sdfj,,x in the spin label (if any), 
will be referred to as magnetic field B. We also present 
some results for the IRLM, for a single channel of spinless 
electrons (A^c = 1): 

Hf^^' = H,, + H,^x+H}^'''' ■ (14) 

In this paper, we focus on the case that Hi and Hp 
differ only in the local level positions (ed^.i 7^ £dfi,F)- It 
is emphasized, however, that our methods are equally 
applicable for differences between initial and final values 
of any other parameters, including the case that the in- 
teractions are channel-specific, e.g. J^fi^j,' t^A'M'^'iA'^rfA'' 

B. AO on Wilson chains 

Wilson discretized the spectrum of Hh on a logarith- 
mic grid of energies zLDA~'^ (with A > 1, /c = 0, 1, 2, . . . ), 
thereby obtaining exponentially high resolution of low- 
energy excitations. He then mapped the impurity model 
onto a semi-infinite "Wilson tight-binding chain" of sites 
fc = to c», with the impurity degrees of freedom coupled 
only to site 0. To this end, he made a basis transforma- 
tion from the set of sea operators {c^^i} to a new set 
{fkp,}, chosen such that they bring H^ into the tridiago- 
nal form 

Afc 00 

^b^EE^^(/l4-i..+h-c-). (15) 

^=1 k=i 
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The hopping matrix elements tk oc DK~^/'^ decrease ex- 
ponentially with site index k along the chain. Because 
of this separation of energy scales, the Hamiltonian can 
be diagonalized iteratively by solving a Wilson chain of 
length k (restricting the sum in Eq. (15) to the first k 
terms) and increasing k one site at a time: starting with a 
short Wilson chain, a new shell of many-body eigenstates 
for a Wilson chain of length fc, say \s)k, is constructed 
from the states of site k and the Mk lowest-lying eigen- 
states of shell fc — 1. The latter are the so called kept 
states \s)^_-^ of shell fc — 1, while the remaining higher- 
lying states \s)^_i from that shell are discarded. 

The typical spacing between the few lowest-lying states 
of shell fc, i.e. the energy scale dE^, is set by the hopping 
matrix element to the previous site, hence 



dEk ~ ife oc i?A-^/2 



(16) 



Now, for a noninteracting Fermi sea with N particles, 
the mean single-particle level spacing at the Fermi energy 
scales as dE oc D/N. This also sets the energy scale for 
the mean level spacing of the few lowest-lying many-body 
excitations of the Fermi sea. Equating this to Eq. (16), 
we conclude that a Wilson chain of length fc represents 
a Fermi sea with an actual size L N , i.e. an effective 
number of electrons TV, that grows exponentially with fc. 



N oc A''^/2 



(17) 



Now consider two impurity models that differ only in 
their local terms Hd,x, and let \Gx)k be the ground 
states of their respective Wilson chains of length fc, ob- 
tained via two separate NRG runs.^ Combining Ander- 
son's prediction (1) and Eq. (17), the ground state over- 
lap is expected to decay exponentially with k as 



|fc(Gi|GF),.| « A- 



with 



A2 



-ak 



log A 



(18) 



(19) 



Thus the AO exponent can be determined by using NRG 
to directly calculate the l.h.s. of Eq. (18) as function of 
chain length fc, and extracting Aao from the exponent a 
characterizing its exponential decay with fc. 

For noninteracting impurity models {U — U' — 0), 
a finite Wilson chain represents a single-particle Hamil- 
tonian for a finite number of degrees of freedom that 
can readily be diagonalized numerically, without the need 
for implementing NRG truncation. The ground state is 
a Slater determinant of those single-particle eigenstates 
that are occupied in the Fermi sea. The overlap {Gi\G-p) 
is then given simply by the determinant of a matrix whose 
elements are overlaps between the I- and F- versions of the 
occupied single-particle states. It is easy to confirm nu- 



merically in this manner that {Gi\Gy) 



, leading 



focus on interacting models henceforth, that require the 
use of NRG. 

In the following three subsections, we discuss several 
technical aspects needed for calculating AO with NRG 
on Wilson chains. 



C. Ground state overlaps 

The calculation of state space overlaps within the NRG 
is straightforward in principle^' especially considering 
its underlying matrix product state structure^^^^^. Now, 
the overlap in Eq. (18) that needs to be calculated in this 
paper, is with respect to ground states as function of Wil- 
son chain length fc. As such, two complications can arise, 
(i) For a given fc, the system can have several degenerate 
ground states {\s)f : s € G}, with the degeneracy dx,k 
typically different for even and odd fc. (ii) The symme- 
try of the ground state space may actually differ with 
alternating fc between certain initial and final configura- 
tions, X e {I, F}, leading to strictly zero overlap there. 
A natural way to deal with (i), is to essentially average 
over the degenerate ground-state spaces, while (ii) can 
be ameliorated by partially extending the ground state 
space to the full kept space, {\s)^ : s £ K}, as will be 
outlined in the following. 

The djf.fe-fold degenerate ground state subspace is de- 
scribed by its projector, written in terms of the fully 
mixed density matrix. 



PG,k 



ax.k 



Ik k\ 



(20) 



It is then convenient to calculate the overlap of the 
ground state space as follows. 



= tr|^fc(PG,/c) 



<h,k 51 $Z l'^^ 



s\s 



vF|2 
Ik 1 



(21) 



to the expected AO in the limit k ^ oo. We will thus 



where tr^j,(-) refers to the trace over the kept space at 
iteration fc of the final system. The final expression can 
be interpreted, up to the prefactor, as the square of the 
Frobenius norm of the overlap matrix ^(sls')!" between 
the NRG states s G G and s' G K at iteration fc for the 
initial and final Hamiltonian, respectively. 

Note that the specific overlap in Eq. (21), as used 
throughout later in this paper, not only includes the 
ground space of the final system at iteration fc, but rather 
includes the full kept space of that system. Yet each such 
overlap scales as e""*^, with the same exponent a for all 
combinations of s and s' , because (i) the states \s)\ with 
s e G are taken from the initial ground state space, and 
(ii) the states \s')J. with s E K from the final kept shell 
differ from a final ground state only by a small number of 
excitations. Therefore, Eq. (21) is essentially equivalent, 
up to an irrelevant prefactor, to strictly taking the over- 
lap of ground state spaces as in ^^^(fc) = trg^(/3g,^). 
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Figure 1: (Color online) Anderson orthogonality for the 
spin-degenerate standard SIAM for a single lead [Eq. (10), 
/i G {tjilli with /^-independent parameters £d and F for Hi 
and Hy as specified in the panel (the full e^-dependence of 
Aao for fixed e\ is analyzed in more detail in Fig. 5) — Sev- 
eral alternative measures for calculating the AO-overlap are 
shown, using Zppi(k) in Eq. (22) with P''^ G {G,K, oo], as 
defined in the text. All overlaps are plotted for even and odd 
iterations separately to account for possible even-odd behav- 
ior within the Wilson chain (thin solid lines with dots, and 
dashed lines, respectively, while heavy symbols identify lines 
with corresponding legends). If even and odd data from the 
same Zppi{k) do not lie on the same smooth line, the com- 
bined data is also plotted (light zigzag lines) as guides to the 
eye. For large k, all AO-overlaps exhibit exponential decay of 
equal strength. Separate fits of e^~"* to even and odd sectors 
are shown as thick solid lines, whose lengths indicate the fit- 
ting range used. The values for A^o extracted from these fits 
using Eq. (19) are in excellent agreement with the displaced 
charge A^j^, as expected from Eq. (8). The relative error is 
less than 1% throughout, with the detailed values specified 
in the legend, and {4Q/lnA) representing the averaged value 
w.r.t. the four measures considered. 

This will be shown in more detail in the following. In 
particular, the overlap in Eq. (21) can be easily general- 
ized to 

z|,p,(A;) = trF,_,(pip,), 0<z|,p,(/c) < 1 (22) 

where P^'^ £ {G,K, oo} represents the ground state 
space, the full kept space, or the ground state taken at 
fc — >■ cx) with respect to either initial or final system, 
respectively. The overlap Zpp,{k) in Eq. (22) then rep- 
resents the fully-mixed density matrix in space P of the 
initial system traced over space P' of the final system, 
all evaluated at iteration k. 

A detailed comparison for several different choices of 
Zpp,{k), including ZQQ{k), is provided in Fig. 1 for the 
standard SIAM with /j, G {t^i})- The topmost line (iden- 
tified with legend by heavy round dot) shows the overlap 
Eq. (2) used as default for calculating the overlap in the 



rest of the paper. This measure is most convenient, as 
it reliably provides data with a smooth fc-dependence for 
large k, insensitive to alternating fc-dependent changes of 
the symmetry sector and degeneracy of the ground state 
sector of Hx.k (note that the exact ground state sym- 
metry is somewhat relative within the NRG framework, 
given an essentially gapless continuum of states of the 
full system). The overlap zgg (data marked by trian- 
gle) gives the overlap of the initial and final ground state 
spaces, but is sensitive to changes in symmetry sector; 
in particular, for fc < 28 it is nonzero for odd iterations 
only. The reason why it can be vanishingly small for cer- 
tain iterations is, in the present case, that the initial and 
final occupancies of the local level differ significantly, as 
seen from the values for {n\^^) and (n^^^) specified in the 
panel. Therefore, initial and final ground states can be 
essentially orthogonal, in the worst case throughout the 
entire NRG run. Nonetheless, the AO-exponent is ex- 
pected to be well-defined and finite, as refiected in zqk- 

The AO-measure zkk (data marked by star) is smooth 
throughout, and although it is not strictly constrained to 
the ground state space at a given iteration, in either ini- 
tial or final system, it gives the correct AO exponent, 
the reason being the underlying energy scale separation 
of the NRG. Finally, Zoo,k = Trf^fc{/5^^} (data marked 
by squares) refers to an AO-measure that calculates the 
overlap of the ground state space of an essentially infinite 
initial system (i.e. k — >■ oo, or in practice, the last site 
of the Wilson chain), with the kept space at iteration 
k of the final system. Since the latter experiences k- 
dependent even-odd differences, whereas the initial den- 
sity matrix ppo^ is independent oik, Zoo,k exhibits rather 
strong /c-dependent oscillations. Nevertheless, their en- 
velopes for even and odd iterations separately decay with 
the same exponent a as the other AO-measures. 

In summary. Fig. 1 demonstrates that all AO- 
measures decay asymptotically as e^~°'^, as expected 
from Eq. (18), with the same exponent a, independent 
of the details of the construction. These details only af- 
fect the constant prefactor A, which is irrelevant for the 
determination of Aao ■ 



D. Channel-specific exponents from chains of 
different lengths 

Equation (6) expresses the exponent Aao of the full 
system in terms of the AO exponents Aao.m '^^ indi- 
vidual channels. This equation is based on the assump- 
tion (whose validity for the models studied here is borne 
out by the results presented below) that for distances 
sufficiently far from the dot, the asymptotic tail of the 
ground state wave function factorizes, in effect, into in- 
dependent products, one for each channel /j,. This can be 
exploited to calculate, in a straightforward fashion, the 
individual exponent Aac/^ for a given channel /i: one 
simply constructs a modified Wilson chain which, in ef- 
fect, is much longer for channel /i than for all others. 
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The overlap decay for large k is then dominated by that 
channel. 

To be explicit, the strategy is as follows. First we need 
to determine when a Wilson chain is "sufhciently long" to 
capture the aforementioned factorization of ground state 
tails. This will be the case beyond that chain length, say 
ko, for which the NRG energy flow diagrams for the kept 
space excitation spectra of the original Hamiltonians Hi 
and Hf are well converged to their T — fixed points 
values. To calculate Aao.^h the AO exponent of channel 
fi, we then add an artificial term to the Hamiltonian that 
in effect depletes the Wilson chain beyond site fcg for all 
other channels ^ /k, by drastically raising the energy 
cost for occupying these sites. This term has the form 



-"art 



(23) 



with C ^ 1. It ensures that occupied sites in the chan- 
nels v fx have much larger energy than the original 
energy scale tk, so that they do not contribute to the 
low-energy states of the Hamiltonian. We then calculate 
a suitable AO-measure (such as zgk) using only fc- values 
in the range k > k^. From the exponential decay found 
in this range, say ^ e^°'^''', the channel-specific AO ex- 
ponent can be extracted, cf. Eq. (19), 



log A 



(24) 



This procedure works remarkably well, as illustrated 
in Fig. 2 for the spin-asymmetric single-lead SIAM of 
Eq. (13) (with Nc = 2, e {t,i})- Indeed, the values 
for Aao,m ^^"^ ^Ao displayed in Fig. 2 fulfill the addition 
rule for squared exponents, Eq. (6), with a relative error 
of less than 1%. 



E. Displaced charge 

The displaced charge Ach,^ defined in Eq. (7) can be 
calculated directly within NRG. However, to properly ac- 
count for the contribution from the Fermi sea, Asoa,/i, a 
technical difficulty has to be overcome: the Hamiltonians 
considered usually obey particle conservation and thus 
every eigenstate of H is an eigenstate of the total number 
operator, with an integer eigenvalue. Consequently, eval- 
uating Eq. (4) over the full Wilson chain always yields an 
integer value for the total Ach,^. This integer, however, 
does not correspond to the charge within the large but 
finite volume l^argc that is evoked in the definition of the 
displaced charge. 

To obtain the latter, we must consider subchains of 
shorter length. Let 
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Figure 2: (Color online) AO exponents for the standard 
spin-degenerate SIAM with spin-asymmetric hybridization 
[Eq. (13), with fi G {tii}] as functions of £d,F (all other 
parameters are fixed as specified in the panel) - The ver- 
tical dashed line indicates ed,i/U — —0.5; at this line the 
initial and final Hamiltonians are identical, hence all expo- 
nents vanish. The squared AO exponents for the individual 
channels, A^o -|- (squares) and A^o i (dots), were calculated 
from Eq. (24). Their sum agrees (with a relative error of less 
than 1%) with A^o calculated from Eq. (19) (down- and up- 
ward pointing triangles coincide), confirming the validity of 
the addition rule for squared exponents in Eq. (6). 



centered on the dot, whose size grows exponentially with 
increasing k. The contribution from channel /i of the 

(k) 

Fermi sea to the displaced charge within ' is 



AL'L^(Gf| 



'^sca,/^ |Gf, 



-(G: 



sea, /I 



|Gi) 



(26) 



where |Gi) and |Gf) are the initial and final ground states 
of the full-length Wilson chain of length N {> k). 

Figure 3 shows Ascl for the spinless IRLM of Eq. (14), 

(k) 

where we dropped the index since A^c = 1- Asca 
exhibits even-odd oscillations between two values, say 
Ag™" and A°^^, but these quickly assume essentially 
constant values over a large intermediate range of k- 
values. Near the very end of the chain they change again 
rather rapidly, in such a way that the total displaced 
charge associated with the full Wilson chain of length A^, 



A 



(JV) 
ch 



AN) 



Adot, is an integer (see Fig. 3), because 
the overall ground state has well-defined particle number. 
Averaging the even-odd oscillations in the intermediate 
regime yields the desired contribution of the Fermi sea 
to the displaced charge, A^ea = 5(A|™" + A"^^). The 
corresponding result for Ach = Agea + Adot is illustrated 
by the black dashed line in Fig. 3. 



sea,/^ 



(25) 



fc'=0 



count the charge from channel /i sitting on sites to fc. 

(k^ 

These sites represent, loosely speaking, a volume V^^.^^^ 



IV. RESULTS 

In this section, we present results for the single channel 
interacting resonant level model [Eq. (14)], and for single- 
lead and two-lead Anderson impurity models [Eq. (13)]. 
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Figure 3: (Color online) Determination of Ach, for the inter- 
acting resonant level model of Eq. (14), for a single specific 
set of parameters for Hi and Hf , specified in the figure legend 
(the Ed.F-dependence of Aao for fixed ed,i is analyzed in more 
detail in Fig. 4. We obtain Ach (dashed line) by calculating 
Adot and averaging the results for even and odd k. To 
reduce the infiuence of chain's boundary regions, we take the 
average over the region between the vertical dashed lines. 



These examples were chosen to illustrate that the vari- 
ous ways of calculating AO exponents by NRG, via Aaoi 
Aph or Ach, are mutually consistent with high accuracy, 
even for rather complex (multi-level, multi-lead) models 
with local interactions. In all cases, the initial and final 
Hamiltonians, Hi and Hp, differ only in the level posi- 
tion: Sd,i is kept fixed, while e^.F is swept over a range of 
values. This implies different initial and final dot occu- 
pations n^^ x = {G x\f>'d^L\G x) 1 and hence different local 
scattering potentials, causing AO. 

AO exponents are obtained as described in the previ- 
ous sections: we calculate the AO-measure zgk{J^) using 
Eq. (2), obtaining exponentially decaying behavior (as in 
Fig. 1). We then extract a by fitting to e""*^ and de- 
termine Aao via Eq. (19). In the figures below, the re- 
sulting A^Q is shown as function of £d/i,F, together with 
A^j^, and also Apj^ in Fig. 4. The initial dot level position 
£d^l.l is indicated by a vertical dashed line. When Ed/^.p 
crosses this line, the initial and final Hamiltonians are 
identical, so that all AO exponents vanish. To illustrate 
how the changes in edp,F affect the dot, we also plot the 
occupancies nd^,F of the dot levels. 



A. Interacting resonant-level model 

We begin with a model for which the contribution of 
the Fermi sea to the displaced charge is rather important, 
namely the spinless fermionic interacting resonant level 
model [Eq. (14), Nc = !]■ The initial and final Hamiltoni- 
ans, ijiRLM g^j^jj ^IRLM^ differ only in the level position: 
the initial one is kept fixed at e^; j = 0, while the final one 
is swept over a range of values, e^^F G [—1, !]• The results 
are shown in Fig. 4. The final dot occupancy nd,F (heavy 
dots) varies from ~ 1 to ~ 0, and Adot — "-d.F — f^d.i 
(dashed line) decreases accordingly, too. The total dis- 
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Figure 4: (Color online) Verification that Aao = Aph = Ach 
[Eq. (5)] for the spinless fermionic interacting resonant level 
model [Eq. 14]. All quantities are plotted as functions of ed,F, 
with all other parameters fixed (as specified in the panel). 
The vertical dashed line indicates ed,i/U' — 0. Heavy dots 
indicate the final occupation of the dot, rid. The exponent 
Aao (light solid line) agrees well with Aph and Ach (trian- 
gles), with relative errors of less than 1%. The local and 
Fermi-sea contributions to the displaced charge Ach are plot- 
ted separately, namely Adot (dashed line) and Asca (dashed 
dotted). The latter is determined according to the procedure 
illustrated, for ed,F/U' = —1.75, in Fig. 3. 



triangles), decreases by a smaller amount, since the de- 
pletion of the dot implies a reduction in the strength of 
the local Coulomb repulsion felt by the Fermi sea, and 
hence an increase in Agca (dash-dotted line) . Throughout 
these changes Aao , ^ph and Ach mutually agree with er- 
rors of less than 1%, confirming that NRG results comply 
with Eq. (5) to high accuracy. 



B. Single-impurity Anderson model 

Next we consider the standard spin-degenerate SIAM 
for a single lead [Eq. (13), ^ G {t,i}] with e^^fi = £d and 
r^t — r. This model exhibits well-known Kondo physics, 
with a strongly correlated many-body ground state. 

In this model, the dot and Fermi sea affect each other 
only by hopping, and there is no direct Coulomb inter- 
action between them {U' = 0). Hence, the contribution 
of the Fermi sea to the displaced charge is nearly zero, 
Agca — 0. Apart from very small even-odd variations 
for the first ~ 35 bath sites corresponding to the Kondo 
scale, the sites of the Wilson chain are half-filled on av- 
erage to a good approximation. Therefore Agca ^ Adot 
(explicit numbers are specified in the figure panels; see 
also Fig. 1), so that Ach,^ in Eq. (7) is dominated by the 
change of dot occupation only,^'' 



Ac\ 



^dot 



(27) 



placed charge, Ach 



idot 



Asca (downward-pointing 



As a consequence, despite the neglect of Agca in some 
previous works involving Anderson impurity models, the 
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Friedel sum rule (Aph = Ach) was nevertheless satisfied 
with rather good accuracy (typically with errors of a few 
%). However, despite being small, Agga in practice is on 
the order of |Asoa| < T/^ and thus finite. Therefore the 
contribution of Agoa to Ach will be included throughout, 
while also indicating the overall smallncss of Agca- In 
general, this clearly improves the accuracy of the consis- 
tency checks in Eq. (5), reducing the relative errors to 
well below 1%. 

The Anderson orthogonality is analyzed for the SIAM 
in detail in Fig. 5. The initial system is kept fixed at the 
particle-hole symmetric point, ed,i = —U/2 (indicated 
also by vertical dashed line in panel a), where the initial 
ground state is a Kondo singlet. The final system is swept 
from double to zero occupancy by varying s^^-p/U from 
—2 to 1. The final ground state is a Kondo singlet in the 
regime nd;j,F — 1/2, corresponding to the intermediate 
shoulder in panel (a). Panel (b) shows the AO-measure 
ZGK{k) as function of k, for a range of different values 
of Ed.F- Each curve exhibits clear exponential decay for 
large k (as in Fig. 1), of the form e^~"^^ . The prefactor, 
parameterized by A, carries little physical significance, 
as it also depends on the specific choice of zpp'; its de- 
pendence on £d.F is shown as thick gray dashed line in 
panel (a), but it will not be discussed any further. In 
contrast, the decay exponent a directly yields the quan- 
tity of physical interest, namely the AO exponent A^q 
via Eq. (19). Panel (a) compares the dependence on e^.F 
of A^Q (dashed line) with that of the displaced charge 
A^jj (light thick line), that was calculated independently 
from Eqs. (7) and (8). As expected from Eq. (5), they 
agree very well: the relative difference between the two 
exponents A^q and A^j^ is clearly below 1% through- 
out the entire parameter sweep, as shown in the inset of 
Fig. 5(b). 

The contribution of the Fermi sea to the displaced 
charge is close to negligible, yet finite throughout (black 
line in panel a). Overall, Agea ^ 0.0037, as indicated 
in Eq. (27). Nevertheless, by including it when calculat- 
ing Ach, the relative error ^A^ is systematically reduced 
from a few percent to well below 1% throughout, thus 
underlining its importance. 

C. Multiple Channels and Population switching 

Figure 6 analyzes AO for lead-asymmetric two-level, 
two-lead SIAM models, with Hamiltonians of the form 
Eq. (13) (explicit model parameters are specified in the 
panels). Panel (a) considers a spinless case (A^c — 2, 
/i = je{l,2}), whose dot levels have mean energy at 
fixed splitting (5, 

Sdi - £d - 5/2, £d2 = £d + 5/2. (28a) 

Panel (b) considers a spinfuU case, [Nc = 4, ^ = (jct) 
with j £ {1, 2}, a S {f, i}), where both the lower and up- 
per levels have an additional (small) spin splitting B <^ 5, 

£djt = £dj + B/2, Edj^ = Edj - B/2 . (28b) 




-2 -1 1 20 40 60 80 

f'U Wilson shell k 



Figure 5: (Color online) Anderson orthogonality for the 
single-lead, spin-symmetric SIAM [Eq. (13), with parameters 
as specified in the legend] - The energy of the d-level of the 
final system £d,F is swept past the Fermi energy of the bath, 
while that of the initial reference system is kept fixed in the 
Kondo regime at £d,i = —U/2, indicated by vertical dashed 
line in panel (a) and in inset to panel (b). Panel (a) shows, as 
function of ed,F, the dot occupation per spin rid^ (dotted solid 
line), the contribution to the displaced charge by the Fermi 
sea, Asea^ (thin black line), the displaced charge A^h (light 
solid line), and the parameters of the large-fc exponential de- 
cay e^""* of ZGK{k) as extracted from panel (b), namely A 
(thick dashed line) and Aao (dark dashed line), derived from 
a via Eq. (19). Panel (b) shows the AO-measure ZGK{k) 
in Eq. (2) (light lines) for the range of Ed.F values used in 
panel (a). The heavy lines shown on top for > 64 are ex- 
ponential fits, the results of which are summarized in panel 
(a). The inset shows the relative error in the AO exponents 
(5A^ = (Aao ~ / ^ti^ the deviation between the light 
solid and dark dashed curve in panel (a); this error is clearly 
less than 1% over the full range of £d analyzed. 

Charge is conserved in each of the Ac channels, since 
these only interact through the interaction on the dot. 
In both models, the upper level 2 is taken to be broader 
than the lower level 1, r2 > Fi (for detailed parameters, 
see figure legends). As a consequence, ■^""^'^ these models 
exhibit population switching: when e^.F is lowered (while 
all other parameters are kept fixed) , the final state occu- 
pancies of upper and lower levels cross, as seen in both 
panels of Fig. 6. 

Consider first the spinless case in Figure 6(a). The 
broader level 2 shows larger occupancy for large posi- 
tive Ed.F- However, once the narrower level 1 drops suffi- 
ciently far below the Fermi energy of the bath as Ed.F is 
lowered, it becomes energetically favorable to fill level 1, 
while the Coulomb interaction will cause the level 2 to be 
emptied. At the switching point, occupations can change 
extremely fast, yet they do so smoothly, as shown in the 
zoom in the inset to panel (a). 

Similar behaviour is seen for the spinfull case in Figure 
6(b), though the filling pattern is more complex, due to 
the nonzero applied finite magnetic field B (parameters 
are listed in the legend). The occupations ridia of the 
narrower level 1 show a strong spin asymmetry, since the 
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Figure 6: (Color online) Anderson orthogonality for a spin- 
less (panel a) and spin-full (panel b) two-lead SIAM, with dot 
levels of unequal width and a split level structure as defined 
in Eq. (28) (all relevant model parameters are specified in the 
legends) - In both cases, the higher level 2 is broader than 
the lower level 1 (r2 > Fi), leading to population switch- 
ing as function of the average final level energy Ed, p. The 
fixed value of ed,i is indicated by the vertical dashed line. 
The inset to panel (a) shows a zoom into the switching re- 
gion, clearly demonstrating that population switching occurs 
smoothly. For panel (b), a finite magnetic field B causes a 
splitting between spin-up and spin-down levels, resulting in 
a more complex switching pattern. In both panels A\q and 
Ach agree very well throughout the sweep, with a relative 
error JA^ well below 1%. 



magnetic field is comparable, in order of magnitude, to 
the level vi^idth {B = ri/2). This asymmetry affects the 
broader level 2, which fills more slowly as is lowered. 
Due to the larger width of level 2, the asymmetry in its 
spin-dependent occupancies is significantly weaker. As 
in panel (a) , population switching between the two levels 
occurs: as the narrower level 1 becomes filled, the broader 
level 2 gets depleted. 

The details of population switching, complicated as 
they are (extremely rapid in panel (a); involving four 



channels in panel (b)), are not main point of Fig. 6. In- 
stead, its central message is that despite the complexity 
of the switching pattern, the relation A^q = A^^ is satis- 
fied with great accuracy throughout the sweep (compare 
light thick and dark dashed lines). Moreover, since Ach 
was calculated by adding the contributions from separate 
channels according to Eq. (8), this also confirms the ad- 
ditive character of AO exponents for separate channels. 

As was the case for the single-channel SIAM discussed 
in Sec. IV B above, a direct interaction between dot and 
Fermi sea is not present in either of the models considered 
here {U' = 0). Consequently, the displaced charge Ach is 
again dominated by Adot, with Agca ^ Adot (c/. Eq. 27). 
Specifically, for the spinless or spinfuU models, we find 
Agca < 0.019 or 0.011, respectively, for the entire sweep. 



V. SUMMARY AND OUTLOOK 

In summary, we have shown that NRG offers a straight- 
forward, systematic and self-contained way for studying 
Anderson orthogonality, and illustrated this for several 
interacting quantum impurity models. The central idea 
of our work is to exploit the fact that NRG allows the 
size-dependence of an impurity model to be studied, in 
the thermodynamic limit of iV — > oo, by simply studying 
the dependence on Wilson chain length k. Three dif- 
ferent ways of calculating AO exponents have been ex- 
plored, using wave-function overlaps (Aao)> changes in 
phase shift at the Fermi surface (Aph), and changes in 
displaced charge (Ach)- The main novelty in this paper 
lies in the first of these, involving a direct calculation of 
the overlap of the initial and final ground states them- 
selves. This offers a straightforward and convenient way 
for extracting the overall exponent Aao- Moreover, if 
desired, it can also be used to calculate the exponents 
Aao./^ associated with individual channels, by construct- 
ing a Wilson chain that is longer for channel /i than for 
the others. We have also refined the calculation of Ach, 
by showing how the contribution Agea of the Fermi sea 
to the displaced charge can be taken into account in a 
systematic fashion. 

The resulting exponents Aao, ^ph and Ach agree ex- 
traordinarily well, with relative errors of less than 1%. 
This can be achieved using a remarkably small number 
of kept states Mk- For example, for the spinfull SIAM 
analyzed below, a better than 5% agreement can be ob- 
tained already for Mk > 32. (For comparison, typically 
Mx = 250 is required to obtain an accurate description 
of the Kondo resonance of the d-level spectral function 
in the local moment regime of this model.) 

Our analysis has been performed on models exhibiting 
Fermi liquid statistics at low temperatures. As an out- 
look, it would be interesting to explore to what extent 
the non-Fermi liquid nature of a model would change 
AO scaling properties, an example being the symmetric 
spinful two-channel Kondo model. 

Finally, we note that non-equilibrium simulations of 
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quantum impurity models in the time-domain in response 
to quantum quenches are a highly interesting topic for 
studying AO physics in the time domain. The tools 
to do so using NRG have become accessible only rather 
recently. One considers a sudden change in some 
local term in the Hamiltonian and studies the subsequent 
time-evolution, characterized, for example, by the quan- 
tity (Gi|e^*'^''*|G'i). Its numerical evaluation requires the 
calculation of overlaps of eigenstates of Hi and Hp. The 
quantity of present interest, [(GiIGf)!, is simply a par- 
ticular example of such an overlap. As a consequence, 
the long-time decay of {Gi\e^^^^^\Gi) is often governed 
by Aaoi too,'^'"" showing power-law decay in time with 
an exponent depending on Aao- This will be elaborated 



in a separate publication. 
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